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Abstract. Given an undirected graph G = (V, E), and a designated vertex q e V, the notion of 
a G-parking function (with respect to q) was independently developed and studied by various 
' authors, and has recently gained renewed attention. This notion generalizes the classical notion 

, of a parking function associated with the complete graph. In this work, we study properties 

}_( ' of maximum G-parking functions and provide a new bijection between them and the set of 

, spanning trees of G with no broken circuit. As a case study, we specialize some of our results to 

the graph corresponding to the discrete «-cube Q„. We present the article in an expository self- 
contained form, since we found the combinatorial aspects of G-parking functions somewhat 
scattered in the literature, typically treated in conjunction with sandpile models and closely 
related chip-firing games. 

o 
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1. Introduction 

m 

> 

'sj- I The classical parking functions provide a bijective correspondence between the spanning 

trees of the complete graph K„ and certain integer- valued functions on the vertices of K„. A 

notion of parking functions corresponding to the spanning trees of an arbitrary graph G is more 

O ' recent and has been independently developed in physics and combinatorics. It was introduced 

00 ' 

1^ I by Bak, Tang and Wisenfeld [jT|| as a self-organized sandpile model on grids, and was general- 

ized to arbitrary graphs by D. Dhar [[TTII . See Definition 12.11 below for the precise definition of 
^ ■ a G-parking function, associated with a connected graph G. 

This notion is already rather powerful; besides generalizing the classical parking function 
from Kn to an arbitrary graph, it has been investigated in the context of chip-firing games 
|[8l |29l [301 and the Tutte polynomial ^ [HI in discrete mathematics, and also investigated in 
algebra and related fields |[2l[T6l|3Tl. However, some of the combinatorial aspects of this topic 
appear somewhat scattered in the literature. 

Several fundamental results concerning the recurrent configurations of chip-firing can be 

derived without the chip-firing context and terminology. For this reason, we shy away from 

introducing and discussing the chip-firing terminology. Instead, in this article we describe 

various interpretations of the G-parking functions in the most elementary combinatorial ways. 
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Using a natural partial order < on the set P(G, q) of parking functions, we consider the maxi- 
mal elements in this poset (!P(G, q), <) . Much of our focus in this paper is on understanding 
the properties of such maximal parking functions. The first result we describe (see Theo- 
rem 14.11) provides a new bijection between the maximal parking functions in the poset and 
the set J?l(G; q) of acyclic orientations of G with a unique source at q. En route, we describe 
what we call an Extended Dhar algorithm (since it is an extension of an algorithm due to Dhar 
IfTTl to recognize G-parking functions) in providing an acyclic orientation corresponding to a 
maximal parking function. We review various combinatorial consequences and algebraic con- 
nections of this correspondence. For example, using known results (namely those of Greene 
and Zaslavsky [{241 and more recent work of Gebhard and Sagan we further identify a 

1-1 correspondence between the set of maximal parking functions and the set of spanning trees 
with no "broken circuits," or equivalently, the set of "safe" spanning trees ; see Section |4] for 
the definitions of these terms. In this paper, we provide a much simpler bijection (compared 
to [EH) between the set of safe trees and the set of acyclic orientations with a unique sink (or 
equivalently, a unique source). Furthermore we generalize this bijection to one between all 
spanning trees and all G-parking functions which preserves the bijection between safe trees 
and maximal G-parking functions. We must remark here that other bijective proofs between 
the set of G-parking functions and the set of spanning trees of G (for arbitrary connected G) 
have been given by Chebikin and Pylyavskyy [[T3]|. However, to our knowledge, the simpler 
bijection we report here, in Theorem 14.21 below, and its generalization given in Theorem 14. 6[ 
are indeed new. 

As an additional contribution, we describe a simple way to generate maximal parking func- 
tions in the Cartesian product graph GiOGi, using maximal functions in the (factor) graphs 
Gi and G2. We then specialize our study to understanding the parking functions in the dis- 
crete n-cube 2„ on 2" vertices. By describing certain special constructions of maximal parking 
functions / on Q„, we obtain a natural description of a set, dom(/), of parking functions - 
those dominated, in the partial order given by ^, by a special maximal parking function /. 
Interestingly enough we shall deduce (see Theorem 15.21) that 
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while it is a well-known fact that 
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Recall that (11.21) corresponds to the total number of spanning trees of 2„ (see equation (5.85) 
in [j321), using the matrix-tree theorem and the explicit knowledge of the corresponding eigen- 
values, to help evaluate the determinantal formula. In light of the fact that finding a bijective 
proof accounting for the number of spanning trees of Qn has been open for several years, we 
hope this is a nontrivial step towards such a proof. 

The paper is organized as follows. In Section [2l we review some preliminaries, including 
Dhar's bum criterion, which determines whether a given function is a parking function. In 
Section |3l we show the bijection between maximum parking functions and acyclic orientations 
with a unique source. In Section 14. 1[ we describe our new and simpler bijection between the 
set of acyclic orientations with a unique sink and the set of safe trees. In Section 15. 1[ we 
describe a construction of maximum parking functions on Cartesian products of graphs. In 
Section \5J\ we focus our study on the n-cube Q„, and provide some explicit constructions 
of maximum parking functions and related bounds. In Section [6l for expository purposes we 
review a bijection between diffuse states (introduced in the context of chip-firing) and acyclic 
orientations of a graph. We conclude with some remarks on research in future directions and a 
few open problems in Section |71 

2. G-Parking Functions and Dhar's Burn Criteria 

In this section we recall the definition of a G-parking function and review Dhar's (burning) 
algorithm that can be used to determine whether an integer- valued function on the vertices of 
G is a G-parking function. 

Definition 2.1. For a connected graph G, a G-parking function relative to vertex q e G is a 
function f : V(G)^Z>^i such that f{q) = -1 and for every non-empty A c V{G)\{q}, there 
exists V e A such that < /(v) < dj(v), where d-^v) is the number of edges e = vw with w t A. 

Remark 2.1. Herein, we have modified the standard definition of a G-parking function some- 
what. The function f is now defined on all ofV{G) instead of restricted to simply V(G)\{q} in 
order to improve the compatibility between G-parking functions and Cartesian product graphs 
such as Qn; due to this change, for f to be a G-parking function, f{q) = -1 necessarily. 

Proposition 2.1. If for a function f : V(G)\{q}^Z>o, for every non-empty connected subgraph 
A c G\{q}, there exists v e V{A) such that f{v) < d-^v), then f is a G-parking function. 

Proof. Assume that, for all connected A c G\{q}, that there exists v € V{A) such that /(v) < 
d-^v). Proceeding by contradiction, suppose that there is some disconnected B Q G\{q} such 
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that /(v) > dg{v) for every v e V{B). Consider then any connected component C of B. Since C 
is connected we have, by the hypothesis of the proposition, that /(v) < d-^{v), for some vertex 
V in C. Thus ^^(v) < d(^{v), implying that there is a vertex u in C\B such that v and u are 
connected by an edge in G; otherwise, either /(v) > d(^{v) or /(v) < dg{v). This contradicts the 
choice of C. ■ 

Throughout we assume that the reference vertex q is fixed, and we always consider parking 
functions with respect to this fixed vertex q, without necessarily bringing explicit reference to 
it. 

A natural question to ask is whether a given integer- valued function on the vertices of G can 
easily be tested for being a G-parking function. In the context of the so-called sandpile models, 
Dhar [fTTl provided an algorithm, which can be interpreted as an efficient algorithm to test if a 
given function is a G-parking function. This was observed in lfT3]| . wherein the algorithm was 
reformulated as follows. Let / : V \ {q} ^ 2^. We assume that f(q) = -1. 

Step 1. Mark any unmarked vertex v which has more marked neighbors than /(v). 

Step 2. Repeat Step 1 until no more vertices can be marked. 

Step 3. Declare / to be a G-parking function if and only if all the vertices have been marked. 

We omit the proof of correctness of the algorithm (as an exercise), which follows in a fairly 
straightforward way from the definition of a parking function. 

On the set P(G, q) of parking functions on G with respect to q, there is a natural partial order 
we may define: 

Definition 2.2. Given two parking functions f,g € V{G,q), we say g < f, if g{v) < /(v), for 
all V e V. The maximal elements in this partial order will be referred to as maximal park- 
ing functions. Finally, a parking function with the largest sum is called a maximum parking 
function. 

For f & r(G,q), let \\f\\ := Zvf(v), and dom(f) = {g £ P(G,q) : g < f}. Furthermore, for 
f,g e f (G, q), let f A g be the function on vertices, defined as f A g(y) := min{/(v), ^(v)}, for 
each V e V(G). 

While the following propositions are perhaps folk-lore, the first part of Proposition 12.21 and 
Proposition [23] appear as Lemmas 7 and 5 in [|9]|. 

Proposition 2.2. (a) Let f e P{G,q), and suppose g : y(G)— >Z>_i such that g(q) = -1 and 
< g(v) < f(v), forv eG with v q. Then g e PiG, q). 
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(b) Iff, g e PiG, q), then f Ag e P(G, q) and dom(/ A ^) = dom(/) n dom(g). 
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Proof. While Part (a) is clear from the definition of a G-parking function, it can also be seen 
using Dhar's algorithm. 

Note that / A g{q) = -1 and so f A g e ViG, q), by Part (a). Now f Ag < f and f A g < g, 
so h e dom(/ A g) implies that h e dom(/) n dom(^). Finally consider, h 6 dom(/) n dom(^). 
Then at each vertex v e V(G), h(v) < f{v) and h{v) < g{v), so h{v) < min{/(v), ^(v)} = f Ag(v), 
proving that h e dom(/ A g). 

■ 

Proposition 2.3. For every connected graph G = (V, E), every f e PiG, q), we have \\f\\ < 
\E\ - \ V\. More over, the equality is always achieved. 

Proof. We may prove this by induction on the number n > 1 of vertices G. The base case 
consisting of V = [q] and no edges is trivially true. For the induction case, let n > 2. Given 
/ e P{G, q), let v„ be the last vertex to be marked by the Dhar algorithm. Then consider 
H := G \ {Vn}, the graph obtained by removing v„ and its incident edges. H is connected, 
since Dhar successfully marks all vertices before v„, and more over, the function fn denoting, 
/ restricted to H, is an //-parking function with respect to q e H. Thus we may apply the 
induction hypothesis to H and fn and complete the proof: 

11/11 = f(Vn) + UhW < J] /(V) < d(vn) - 1 + {\E\ - J(v„)) - (\V\ - 1) < \E\ - \V\ , 

where we also used the fact that /(v) < d{v) - 1, for every v and parking function /. 

The proof also suggests that by assigning the maximum possible value, at each step in Dhar's 
marking algorithm, one easily obtains a (maximum) parking function which achieves the upper 
bound. ■ 

Note that the quantity g(G) := l^l - \V\ + 1 is sometimes referred to as the cyclomatic number 
or the Betty number of the graph and due to our convention of assigning f{q) := -1, we have 
g(G) - 1, as the bound in the above proposition. 

Proposition 12.21 also gives us a simple (albeit not necessarily efficient), inclusion-exclusion 
method to relate the set of maximum parking functions with the set of all parking functions. 
But this has to wait until the next section, where we observe another basic fact concerning the 
maximum parking functions. 
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3. Maximal G-Parking Functions and Acyclic Orientations with a Unique Source 

Given a graph G, the notion of an acycUc orientation of the edges of G is classical, with an 
extensive literature. The notion of an acyclic orientation with a unique source at a fixed vertex 
is less well-studied. Let AO(G) denote the set of acyclic orientations of the graph G and let 
J?l(G, q) be the set of acyclic orientations of G with a unique source at vertex q. Finally, let 
MP(G, q) denote the set of maximum G-parking functions. 

Theorem 3.1. There exists a bijection between ^(G,q) and the set MP(G,q) of maximum 
G-parking functions. 

Proof. Given an acyclic orientation O e ^(G, q) with a unique source at q, define the function 
/ = f(0) on the vertices of G: let f(v) be the indegree (in O) of v minus 1. We will show that 
this correspondence provides the necessary bijection. 

In any orientation, the sum of the indegrees equals the number of edges. Hence 2v /(^) = 
l^l - By using Dhar's algorithm, we may show that / is in fact a parking function: starting 
with q, we may repeatedly mark and remove the current set of source(s) in the acyclic orien- 
tation of the remaining graph; since a vertex v with value f(v) becomes a source only when 
all its f(v) + 1 in-neighbors have been marked and removed, the Dhar criterion is satisfied. 
Also observe that the procedure stops only after marking all the vertices, since every acyclic 
orientation has at least one source. Thus f(0) G MP(G, q). 

To see that f{0\) fiOj), whenever Oi 4^ O2, simply recall that an acyclic orientation is 
uniquely determined by its outdegree sequence: starting with the sinks, orient all edges into 
the sinks, remove the sinks, and repeat the process by subtracting one from the outdegrees of 
the neighbors of the sinks. 

The proof will be complete once we establish the onto property, that every maximum parking 
function can be obtained this way. Given a maximum parking function / e MP(G, q), we will 
construct an orientation 0(f) using the following modification of Dhar's algorithm, and will 
show that 0(f) e Jl(G; q), thus essentially providing an inverse map to the above construction. 

The Extended Dhar Algorithm. 

Input: A maximal parking function / e MP(G, q) 

Output: An acyclic orientation 0(f) with a unique source at q. 

Step 1. Start with v = q. Orient all edges out of q. 
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Step 2. If there exists a vertex v which accrued indegree(v) equal to /(v) + 1, mark v and 
orient the remaining edges incident at v outward from v. 

Step 3. Repeat Step 2, until all vertices are marked and all edges are oriented. 

The correctness of the original Dhar algorithm guarantees that all vertices will eventually 
be marked - indeed, the indegree(v) equals the number of neighbors marked before v; thus all 
edges will be oriented, meaning that 0{f) is an orientation of the edges of G. Observe that the 
indegree of a vertex v equals /(v) + 1. Since q is unique with f{q) := -1, it must be that q 
is the unique source. It is also easy to see that 0{f) is acyclic - if there were to be a cycle, 
considering the first vertex in the cycle which was marked, we obtain a contradiction to the 
way the edges were oriented (in Step 2 above) from a marked vertex. I 

Remark 3.1. Upon completion of this work, we discovered (thanks to Matt Baker), that The- 
orem I3.il can also be derived using chip-firing games: As described in ll23l . the notion of 
a so-called diff"use state (introduced by ll25l and see Definition \6.1\ below) helps relate chip- 
firing configurations to acyclic orientations. Also thanks to an anonymous referee of an earlier 
version of this work [|6l, we learned that Theorem \3.1\ is Lemma 10 (under the name o/ allow- 
able orientations) of Biggs where it is mentioned that this in fact goes back to an even 
earlier result of Greene and Zaslavsky Il24l . 

It is now easy to observe the following fact (which appears as Lemma 8 in (HI). 
Corollary 3.2. Every maximal parking function is a maximum parking function. 

Proof. This follows from the proof of correctness of the Extended Dhar algorithm described 
above - if / were maximal, but not maximum, then there must be a vertex in Dhar's marking 
whose indegree is at least /(v) + 2. But then we can increase /(v) by one, and obtain a valid 
parking function, contradicting the maximality of /. I 

We now return to prove the simple result that was promised at the end of the previous section. 

Corollary 3.3. Let G be a finite graph with k maximum G-parking functions. Then there exist 
G -parking functions fx,. . . ,fn such that the number of G-parking functions is 

n 

2+|dom(/;)|, 

!=1 

where n = 2^ - 1 and the sign + or - is uniquely determined by Proposition \2.2\ 
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Proof. Let gi,g2,---,8k e 'P(G,q) be all the maximum G-parking functions. Trivially, every 
non-maximal parking function is dominated by some maximal parking function, and Corol- 
lary |X2l lets us observe that. 

Now we may simply use the inclusion-exclusion formula to count the size of the union of the 
above k sets: 

\'P(G,q)\ = |utidomfe)l 

= |dom(^,)| - Yj Idomfe) n dom(^y)| + ■ • • 

I \<i<j<k 

+(-l/|dom(^i) n dom(g2) n • • • n dom(^^)| 
= ^Idomfe)!- Yj |dom(^; Ag,)| + --- + (-l/|dom(^i Ag2 a--- A^^)|, 

i \<i<j<k 

which, upon using Proposition [2]2l completes the proof. ■ 

It would indeed be interesting to see if the above corollary can be used in making progress 
towards obtaining a bijective proof for the number of spanning trees of Q„, the n-dimensional 
hypercube. In Section |5^ we take a modest step towards it. 

Theorem 3.4. Let G be a simple, connected graph. Then, for a fixed choice of q, G has a 
unique maximum G-parking function if and only ifG is a tree. 

Proof. If G is a tree, then there is only one parking function with respect to any q since G has 
no cycles and, thus, each vertex can have at most one marked neighbor in the Dhar algorithm. 
(Note that this is in fact tautological if one uses the bijection between the parking functions 
and the spanning trees of G.) Hence there is only one maximum G-parking function. 

The other direction is less obvious. However, observe that in light of the bijection established 
in Theorem 13. 1[ it suffices to show the following. Whenever G is connected and contains a 
cycle, then there are at least two acyclic orientations for G, with q as the unique source. This 
is easy to establish (for example, by considering the standard directed acyclic graph (DAG) 
representation of the graph), and we leave the proof as a simple exercise. ■ 

See Remark |4TT] below for another short (but indirect) proof of the above theorem. 



Corollary 3.5. A simple, connected graph G has a unique maximum parking function f if and 
only if the range of f is a subset o/{-l, 0}. 
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Proof. This follows from the observation that any tree has a unique parking function and, for 
any vertex v q, /(v) = 0. ■ 

Corollary 3.6. For every G, WiG, q)\ = 1 if and only if \MF(G, q)| = 1. 



4. BijECTiONS OF Maximal G-parking Functions 

We now augment our Theorem 13.11 above to show a bijection between maximal G-parking 
functions and spanning trees with no broken circuits. Before stating the theorem, we need to 
recall several definitions from the literature. 

The Tutte (or Tutte-Whitney) polynomial of a graph G = (V, E) is the two-variable polyno- 
mial defined as 

(4.1) Tcix^y) = _^(jc - lf(^)-<^) (y - , 

ACE 

where n = |y| and k(A) denotes the number of connected components of the graph on V 
using edgeset A. For A e the chromatic polynomial ^g('^) of ^ graph G is defined as the 
number of proper vertex colorings of G using A colors. For a general variable A, the following 
relation between the chromatic polynomial and the Tutte polynomial is well-known (see e.g. 
mill 1331 121 EH): 

(4.2) xciA)) = (-l)"-(^) /l'^(G) Toil -A,0), 

where n = |y| is the number of vertices and /c(G) is the number of connected components of G. 

Definition 4.1. Given a graph G = iV,E) and an ordering of all the edges of G, a broken 
circuit B Q E is any cycle (of edges) ofG minus the largest (according to the ordering) edge in 
the cycle. 

Note that since every cycle contains (or gives rise to) a broken circuit, a collection of edges 
not containing a broken circuit must necessarily be acyclic. Inspired by the terminology of 
Kenyon and Winkler [|27]| . we call a spanning tree T safe, if it contains no broken circuits. That 
is, for all edges e not in the tree, there is an edge in the unique cycle formed when e is added 
to the tree, which is larger than e. The notion of a broken circuit is more general, and in fact 
explains the classical terminology: in the context of a matroid, an independent set of elements 
of the matroid obtained from a circuit, by removing the largest element (once again, according 
to some a priori global ordering of all the elements) of the circuit. 
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Theorem 4.1. For every undirected, connected graph G, the following quantities are all the 
same. 

• (a) The number of maximum G-parking functions with respect to q. 

• (b) The number of acyclic orientations with a unique source at q. 

• (c) The number of spanning trees with no broken circuits, or equivalently, with zero 
external activity. 

• (d) The coefficient (up to sign) of the A-term in the chromatic polynomial xg{X). 

• (e) The value (up to sign) 7~g(1, 0) of the Tutte polynomial Tcix^y), evaluated at x = I 
and y = 0. 

Proof. Theorem 13711 shows the equivalence of (a) and (b). The equivalence between (c) and (d) 
is classical and is part of Whitney's Broken Circuit theorem Il35l : that the chromatic polynomial 
on n vertices is given by 

n-l 

Xg(^) = Yji-iymrA'"-'-, 

where mj is the number of r-subsets of edges of G which contain no broken circuit. The term 
m,j_i corresponds to (the absolute value of) the coefficient of A; note that the n - 1-subsets under 
consideration being necessarily acyclic, correspond to spanning trees which do not contain a 
broken circuit. 

The equivalence between (b) and (d) is due to Greene and Zaslavsky [|24J . A (direct) bijective 
proof of the equivalence between (b) and (c) is given by Gebhard and Sagan [|2T|. using a 
modification of an algorithm of Blass and Sagan [[T2|. In Section 14.11 we provide a much 
shorter proof of the equivalence between (b) and (c). 

The equivalence between (d) and (e) follows from (|4.2I) . and using k(G) = 1, for a connected 
G. The equivalence between (a) and (e) follows from results of [1171 and [|29l . which confirmed 
a conjecture of Biggs f8| in the context of chip-firing. An inductive proof (using edge deletions 
and contractions) without involving chip-firing is due to Plautz and Calderer [|30l . As described 
in [|30l . the work of Dhar and Lopez provides the following result: 

fenG,q) 

where w(f) = [^(G)! - |y(G)| -I- 1|/||, hence the equivalence of (a) and (e). The results in [29] and 
[[T4| also establish the equivalence between (c) and (e), with the minor modification that broken 
circuits are equivalently described using external acitivities - each broken circuit contributes an 
external activity of one to a spanning tree. ■ 



G-PARKING FUNCTIONS, ACYCLIC ORIENTATIONS AND SPANNING TREES 11 

In addition to the above, Cori and Le Borgne lfT4l describe certain decreasing traversals of 
vertices and edges and a notion of strong edges to provide a bijection between recurrent chip- 
firing configurations (with a fixed "level") and spanning trees with a fixed "external activity". 
While the level corresponds to the sum of the values of a parking function (up to an additive 
shift), the external activity reflects the number of broken circuits, and we refer the interested 
reader to their paper for additional information. 

Remark 4.1. Observe that the nontrivial part of Theorem \3.4\f ollows easily using the equiva- 
lence between (a) and (c). Indeed, let G be a connected graph which contains a cycle of length 
3 or more. Then given a spanning tree which contains no broken circuit, we include an edge 
not in the tree to form a cycle C. Since C is of length 3 or more, there must be an edge, which 
is not the largest edge in the cycle, that can be removed, giving another spanning tree with 
no broken circuits. Hence a ( connected) graph containing a cycle has more than one parking 
function. 

4. 1 . Bijection between Acyclic Orientations with Unique Sink and Safe Trees. In this sec- 
tion we give a shorter proof, of equivalence of (b) and (c) of Theorem 14.11 than the ones 
reported in Gebhard and Sagan [|2T|| and Gioan and Las Vergnas 1221 . Note that the bijection 
in [1211 is not activity -preserving while the one in Il22l is; also see [|28ll . 

For ease of presentation, we will consider orientations with a unique sink rather than a unique 
source - clearly this is equivalent. 

Let cr be any total ordering of the edges of G. Given two edges e and /, we say e is larger 
than / if cr(e) > cr{f). Similarly we say e is smaller than /, if the inequality is otherwise. 
Recall, a spanning tree T of G is called safe with respect to cr if for any edge e iT , there exists 
at least one edge / in the unique cycle inT + e such that / is larger than e. Let T(G, cr) be the 
set of safe trees with respect to cr, and let ^(G, q) be the set of acyclic orientations of G with q 
being the unique sink. 

Theorem 4.2. For any total order cr, there exists a bijection jXa- '■ TiG, cr) ^(G, q) . 

Henceforth, we fix cr and do not write it as a subscript. Before we proceed, we make a few 
more definitions and observations. An arborescence with root ^ is a directed spanning tree 
with all vertices except the root having out-degree exactly 1 , and the root having out-degree 
0. Any spanning tree corresponds to a unique arborescence with root q and henceforth we will 
use the terms interchangeably. Given a vertex /, we let denote the unique directed path from 
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/ to q. Given vertex / and j, we let meet(i, j) be the first vertex in the intersection of P, and 
Pj. That is, the path from / to meet{i, j) and the path from j to meet(i, j) are disjoint except at 
meet(i, j). If j hes on P,, we let j be meet(i, j). Observe that for three vertices /, j, k, either (a) 
meet(i, k) = meet(j, k), or (b) meet{i, k) lies on Pi and meet(j, k) = meet(i, j), or (c) meet(j, k) 
lies on Pj and meet{i, k) = meet{i, j). Given two vertices / and j, we will denote the largest 
edge in the path from / to meet{i, j) as eij and the largest edge from j to meet{i, j) as 6^. If 
j = meet(i, j), we let ep be the null edge. 

We (abuse notation and) say i >a- j if £,7 is larger than ej;. We will also define e,-,- to be a null 
edge. Note that if i >,y j and j >o- k, then going over the three possibilities of meet(i, k) we see 
that / >(;- k. Thus >o- is transitive and induces a total ordering of vertices. We say i dominates 
j, if i >a- j. It is instructive to note that / dominates all vertices in P,. 

We now describe a mapping ju from all arborescences with root q to acyclic orientations 
with unique sink q. We will prove that distinct safe trees lead to distinct arborescences. This 
proves that the mapping fj. restricted to safe trees is one-to-one. Furthermore, given an acyclic 
orientation, we describe a procedure n which takes an acyclic orientation and returns a safe 
tree. Moreover, for any orientation O e J?l(G, q), we have fx{n{0)) = O. This shows that [i is 
onto and thus it is a bijection. 

Arborescence to Acyclic Orientation (ju): Given an arborescence, orient an edge as / 
to 7 if / >o- 7, or vice-versa. By the transitivity of >„- it is clear that the orientation is acyclic. 
Also every vertex dominates the root which therefore is the unique sink. We will call an ori- 
entation so obtained as one induced by the arborescence. The following lemma will show that 
two safe arborescences cannot lead to the same orientation. 

Lemma 4.3. Let T\ and T2 be two distinct safe arborescences and Oi = //(Ti) and O2 = /^(T2). 
Then Oi ^ O2. 

Proof. We prove the contrapositive: suppose Oi = O2 = O, then we show that Ti = T2. 
Consider the trees rooted at q (note that q is the unique sink of O) with edges directed towards 
q. We now show that for each vertex of the graph: 

(*) The unique out-neighbor in Ti is the same as that in T2 which will imply that both trees are 
the same. 

Since O is acyclic with a unique sink, the vertices V can be decomposed as V = (q =: S U 
Si U ... U Sr) for some r > 1, where 5, is the set of vertices which are sinks in the digraph 
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G \ {j'iJoSf. Let Si be the first set (with the least /) to contain a vertex violating (*). Let this 
vertex be denoted / (abusing notation). Let (i, f) and (/, k) be the unique out-neighbors of / 
in Ti and T2 respectively, with j k, j e S j and k e (again abusing notation). Observe 
that j, k < i (that is, S j,Si^ precede Si), by the nature of the decomposition, and the fact that 
= n(Ti) and = n(T2). 

Now consider the undirected cycle using the edges (/, j), (i, k) and the unique paths, Pj from 
j to q and P/, from k to q. Observe that both Pj and P^ are contained in Ti and T2, by the choice 
of i. Also note that the largest edge in this cycle must be either on Pj or on P^, since both Ti 
and T2 are safe! We now get a contradiction - if the largest edge is on Pj, then our definition 
of ju demands that (i, j) be oriented from j to i in Ti ; similarly if it is on Pk, then the edge (/, k) 
be oriented from k to i mT2. □ 

Acyclic Orientations to Safe Trees(;r): Now we describe a procedure to get a safe tree from 
an acyclic orientation O with single sink q. Let d(v) denote the out-degree of vertex v. Note 
that d(q) = 0. For an edge oriented / to j, we say j is a out-neighbor of / and / is an in-neighbor 
of j. At each step we maintain a set of labeled vertices X and an arborescence T spanning X. 
We maintain the invariant that there is no edge from a vertex in Z to a vertex in F := V \X, 
(think of Z as a large sink). Initially, X = {q}. Note that since the orientation is acyclic, at any 
step there is at least one vertex w in F which is a sink in the induced graph G[Y], that is, u has 
no out-neighbors in Y. We pick one such u arbitrarily. Also, since there is a unique sink q, 
this vertex u must have at least one out-neighbor in X. We add w to Z and we connect u to the 
arborescence T as follows. 

Let Xu Q Xhe the set of out-neighbors of u. Let ;c 6 Z„ be the vertex which dominates all 
other vertices in X„ with respect to the current arborescence T. Let Q Z„ be the subset of 
all vertices v such that (u, v) >o- e^v That is, the edge (u, v) is larger than the largest edge in the 
path from x to meet{x, v). Note that x e Wu and thus Wu is non-empty. Connect u to the vertex 
V with the largest (m, v) among all v in Wu- We end when X = V with an arborescence T. 

Lemma 4.4. The arborescence T obtained at the end is a safe tree. Moreover, j^(n(0)) = O, 
forO&Ji(G,q). 

Proof. Consider an edge (j, i) not in the tree. Note that at each step exactly one node is added 
to the arborescence. Also note that for an arc oriented (j, i) in O, i is added before j. Let Xj 
be the set of labeled vertices in the step when j is added to the arborescence. Note that / G Xj. 
Suppose X e Xj was the dominator of Xj and {j, k) was the edge added at this step. 
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Observe that whenever the procedure adds a new vertex j, the vertex j dominates all other 
vertices in the arborescence. This is because {j,k) >a- e^k, that is, the largest edge in the path 
X to meet(x, k) = meet{x, j) and thus j dominates x and so every other vertex. Thus the tree T 
induces the same orientation O. It remains to show that T is safe, that is, (j, i) is not the largest 
edge in the cycle T + (j, i). 

If i e Wj, then by choice of k (based on how the edge (j, k) was added above), (j, k) >a- (j, i) 
and thus (j, i) is not the largest cycle in T + (j, 0- 

If / i Wj, this means (j, /) < Cxi- That is, (j, /) is smaller than the largest edge in the path from 
X to meet(x, i). Also, by definition, (j, k) > e^j that is, (j, k) is larger than the largest edge in the 
path from x to meet{x, k). Let the path from x to meet(x, i) be Qi and that from x to meet{x, k) 
be Qk. 

Now consider meet(i, k). If meet{i, k) does not lie on P^, then meet(x, i) = meet{x, k). Thus, 
(j, k) >cr (j, i). If meet{i, k) lies on P^, then it must be meet(x, i) or meet{x, k). If the latter, then 
Qi ^ Qk, and therefore (j, ^) >o- (7, 0- If meet{i, k) = meet(x, i), then Qk c g,. Either the largest 
edge in 2, lies in Qk and we are done as before; Or, the largest lies in Qj \ Qk, which lies in the 
cycle formed inT + (j, i). Thus, (j, i) is smaller than the largest edge in the cycle ofT + (j, i), 
completing the proof. □ 

Remark 4.2. Thanks to an anonymous referee, we learnt that Biggs and Winkler ifTOl had 
actually given a simple bijection somewhat similar to the one we describe above, but unlike 
ours, their bijection does not preserve the so-called external activity. This is an important 
distinction for us, since we extend the above notion in the next subsection to provide a more 
general bijection. 

4.2. Extension to a bijection of G-Parking Functions and Spanning Trees. In this section 
we generalize the above to a bijection between G-parking functions and all spanning trees of 
G, in such a way that it preserves the bijection between maximal parking functions and safe 
trees. We use the definitions of the previous section and make a few more definitions and claim 
below before demonstrating the bijection. 

Given an arborescence T spanning only a subset of vertices X Q V, and a vertex u i X, 
we make a few definitions and observations which will be useful in our bijection. Let r(M) 
be the neighbors of u and let X„ := Y{u) n X. We now describe an order on the vertices X„ 
(which could be different from the total order of the previous paragraph, but is related), which 
we call the power order of Intuitively, given two neighbors v and w of m in X„, if v is more 
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powerful than w, then in the tree T + {u, v), u would dominate w and the edge {u, w) doesn't 
form a broken circuit with T + {u, v). Moreover, in the tree T + (u, w), either v dominates u or 
the edge (u, v) forms a broken circuit with T + (u,w). 

Let X be the vertex in X„ which dominates all other vertices in X„ with respect to the current 
arborescence T. Let Wu Q be the subset of all vertices v such that (u, v) >o- e„. That is, the 
edge (w, v) is larger than the largest edge in the path from x to meet{x, v). Note that x 6 and 
thus Wu is non-empty. Let v be such that (w, v) is largest among all v in W„. Call v the most 
powerful element of X„. Delete v from X„ and repeat till one gets an order on all vertices of X„. 
We call this order the power order w.r.tT . In Figured] we give an illustrative example. 

q 




Figure 1 . The numbers on the edges correspond to the cr value. The dark edges 
form the tree T and u ^ T. Note that a >a- b >„■ d c mT . However, the 
power order of X„ = {a,b,c,d\ is as follows. Firstly, x = a and Wu = {a,b}. 
Since (u, b) is bigger, b is the most powerful. Second powerful is a. After that 
X = d and Wu = {c, d), and thus the complete power order is {b, a, c, d). 

Claim 4.5. Given a tree T and a vertex u ^ T with neighbors X„ in T, and two vertices v and w 
in Xu- Ifv is more powerful than w, then u dominates w in the tree T + (u, v) and (u, w) doesn 't 
form a broken circuit with T + (u, v). In the tree T + (u, w), either v dominates u or {u, v) forms 
a broken circuit with T + (u, w). 

Proof Consider the time when the power order of v is determined. Let x be the dominator at 
that stage and let Xu the set of neighbors of u remaining. Note w e Xu. Either v e Wu or v = x 
and Wu is empty. If the former, (m, v) is larger than Cxx, and thus inT + (u, v), u dominates x 
which dominates w. Moreover (u, w) is either smaller than (u, v) if w e Wu, or smaller than e^w 
which is smaller than euw since u dominates x. In any case, (u, w) doesn't form a broken circuit 
in r + (m, v). Also, in r + (w, w), {u, v) forms a broken circuit since it is larger than the largest 
edge from v to w. 
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If the latter, then u dominates w inT + {u, v) since v = x dominates w in T . Also since Wu 
is empty, the edge {u, w) is smaller than ey„ implying (m, w) doesn't form a broken circuit in 
T + (a, v). Also, in r + (a, w), v dominates u since {u, w) is smaller than e^w- ^ 

Now we are ready to present the next theorem which gives the desired bijection between 
G-parking functions and spanning trees of G. Let 7~(G) be the set of all spanning trees of G. 
Recall that 7~(G, cr) was the set of safe spanning trees with respect to cr. 

Theorem 4.6. There exists a bijection n : f (G, q) — > 7~(G) such that for all f e V{G, q) which 
is maximal, n{f) is a safe tree. 

As in the proof of Theorem 14.11 we describe mappings ji : 7~(G) — > ^(G, q) and n : 
P(G, q) —> T{G) and show that for any parking function /, jj.{n(f)) = f (Lemma 14.101) and 
for any spanning tree T, n{p.{T)) = T (Lemma [4.1 Ih . Furthermore we show that for a tree T, 
IJ.{T) is a maximal parking function iff T is safe (Proposition l4.7l) . This completes the proof of 
Theorem 14.61 

Spanning Trees to Parking Functions ip): Given a spanning tree T , let E' be the subset of 
edges not in T which form a broken circuit with T . Delete E' from E. Direct the edges of T 
with q as the unique source - that is, all vertices except q have an out-degree of 1 and q has 
an out-degree of 0. Given this arborescence, for any undirected edge (/, j) in £ \ £", orient it 
as (/, j) if / >o- J, or (j, /) if j >a- i. This gives di-graph D. The parking function / := iu{T) is 
defined as /(v) = out-degree£,(v) - 1 in the orientation of the edges of E \ E'. 

Proposition 4.7. / := fx(T) is a maximum parking function ijfT is a safe tree. 

Proof. Note that if T is safe, E' is empty and therefore all the edges are oriented. Thus 
Zv/(^) = 1^1 ~ 1^1 ^rid by Proposition 12.31 / is a maximal parking function. On the other 
hand if T is not safe, then E' is not empty and thus Xv/(^) < 1^1 ~ 1^1 implying / is not a 
maximal parking function. □ 

Parking Functions to Spanning Trees {n)'. Given a valid parking function, we use a modifica- 
tion of Dhar's algorithm to obtain the spanning tree. We maintain a set of vertices X, connected 
via an arborescence T,, with Xq initialized to [q] and Tq is 0. At each step we add one vertex to 
Xi and one edge to T,. In the end we get X„^\ = V and T„^i is the spanning tree returned. We 
describe the i + 1th step. Let r(v) denote the neighbors of v. 
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(1) Let Si^, ■.= {veV\ X,- : |r(v) n X,| > /(v)}. 

(2) For every u e Si+u X„ := r{u) n X,- and let M(u) be the QXJ - f(u))th vertex in the 
power order of X„ with respect to T. Let y„ be the vertices in X„ more powerful than 
M(u) and Eu be the set of - /(w) - 1 edges of the form (u, v) where v 6 y„. 

(3) Note that adding all edges of the form (u, M{u)) to T, gives a new tree T' . In T', let u 
be the vertex in S ,+i which is dominated by all other vertices in S ,+i with respect to T'. 
Add u to get and the edge {u, M{u)) to get r,+i. 

Claim 4.8. Given a parking function f, let the vertices be added in order {q = uq,u\, - • • , m„-i}. 
That is Xi := {uq, • • • , m,}. Then with respect to the tree T,, Uj dominates which dominates 
M;_2 and so on. 

Proof. The proof is by induction. At stage /, let 5, be the set of vertices as defined above. 
Suppose Uj is added at this stage. We will be done if we show w, dominates 

Two cases arise: If is not a neighbor of M;, then m,- must have been in S as well implying 
Ui dominated m,_i in T'._^ and since A/(m,) doesn't change, dominates in T, as well. 

If Ui-i is a neighbor of m,, then either | = /(m,) + 1 and m, connects to the most powerful 
of its neighbors in X„. and thus dominates which is in Or, | > /(w;) + 1, which once 
again implies m, was in 5;-i and moreover, either A/(m,) becomes or still remains more 
powerful than In any case, dominates in T,. □ 

Claim 4.9. For a«j vertex u, the set of edges are precisely the set of edges from u to X„ 
which form broken circuits with T. 

Proof. Consider vertex u connecting to vertex w in X„. From the first part of Claim |43] we see 
that the edges not in do not form broken circuits. Let {u, v) be an edge in Since v is more 
powerful than w, from the second part of Claim |43] we see either (m, v) forms a broken circuit 
with r or V dominates u 'mT . The second possibility is precluded by Claim 1481 □ 

Lemma 4.10. For any parking function f, we have fi{n{f)) = f. 

Proof. Let T be the arborescence formed by rooting the tree nif) at q. Fix a vertex u. From 
Claim 1491 we have that /j. will first remove all the edges in £„. After the removal of these edges, 
u will dominate the remaining f{u) + 1 vertices in Z„ and thus its out-degree will be that. Thus, 
iu(n(f))iu) = out-degree^(M) - 1 = f(u). □ 



Lemma 4.11. For any spanning tree Z, we have n(ju(Z)) = Z. 
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Proof. Let / := //(Z). Abuse notation and call the arborescence obtained by rooting Z at q, 
also Z. Let E' be the set of edges which form broken circuits with Z. Let D be the di-graph 
obtained by orienting the edges of E\E' with respect to Z. Given a subset of vertices X, let 
Z\X\ be the induced sub-forest of X. 

The proof proceeds by induction on the stages of the algorithm computing n. We assume 
at stage /, the current tree of the algorithm, T , is a subtree of the tree Z. That is T = Z[Xj]. 
We also assume that for every vertex in Z,, all its out- neighbors in D are also in X,. These are 
vacuously true at stage 0. We now show that at stage / -i- 1 a vertex connects to tree T using an 
edge of Z and all its out-neighbors of D are in X,. 

Consider vertices mV\ Xj. At least one of these vertices v must have all its f(v) + 1 out- 
neighbors of D in Xi for otherwise we would get a cycle in D. Call this set of vertices S . Note 
that Z[Xi U 5 ] is connected as the vertices in X can only connect to Z using some edge of D. 
Let u* e S he the vertex which is dominated by all other vertices in S in Z[X, U 5 ] . Let {u*, w*) 
be the edge in Z with w* e X, . We claim that the algorithm which computes n also picks u* in 
this stage and w* is M(u). 

Note that S c since each vertex in S has at least f{u) + 1 neighbors in Xj. In fact, we 
show for every vertex u € S , the corresponding w in X, where {u,w) e Z is in fact M{u). Call the 
set of these f{u) + 1 out-neighbors of u in X,, Y^. Observe that any edge of the form (u, v) with 

V i Yu must form a broken cycle with Z. This is because these edges are not directed towards 

V and cannot be directed towards u by the induction hypothesis. Since they form broken cycles 
with Z, these {XJ -{YJ vertices must be more powerful than w with respect to the arborescence 
Z restricted to vertices of X,, that is T by the induction hypothesis. This follows from the 
definition of power. Moreover, w must be powerful than all other vertices of 7,, with respect to 
T since u dominates all these vertices. Thus w is the (IXjJ - \YJ + l)th powerful vertex in X„ 
with respect to T, that is, w = M(u). 

We will be done if we show any vertex u' e S ,+i \ Q dominates some vertex in S with respect 
to Z. If this is the case, then the algorithm would choose the vertex which is dominated by all 
vertices in and it has to be the vertex u*. But this is true since u' has some out-neighbor 
of D in y \ X, - a path following argument shows we must reach a vertex v e 5 from u' using 
edges in D. In other words, u' dominates v with respect to Z. □ 



In Figure [21 below, we give an example of the bijection on a simple 4- vertex 5 -edge graph. 
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Figure 2. We illustrate our bijection on the graph on the top left comer. The 
numbers on the edges is the ordering cr. The graph has 8 spanning trees and they 
are shown beside the graph. The numbers on each spanning tree correspond to 
the 8 possible parking functions of the graph. It is instructive to note that the 
second row of spanning trees correspond to safe trees w.r.t to cr and correspond 
to maximum parking functions - note the sum of the numbers is 1 for the bottom 
4 trees. We now take the first tree and show how it leads to the parking function 
via the function /i and then how vice-versa is obtained via the function n. jj. 
first recognizes the edges which form broken circuits with T - the dotted edge 
is the only one in this case. Then it orients the tree edges towards the root a. It 
orients the edge (c, d) towards d because c >o- d w.r.t T . After the orientation, 
the parking function is found by subtracting 1 from the out-degrees. Going from 
the parking function to the tree, the dark vertices denote the set X, at each step, 
while the dotted lines enclose the set S,. When 5, has size more than 1 (when 
it contains vertices {c,d}), it finds A/(c) = b and M{d) = b according to the 
power-order. It then chooses (d, b) since c dominates d in the tree with both 
(c, b) and {d, b) added. 

We end this section by using the relation between parking function and orientations to prove 
a property about the poset formed by parking functions. Recall given two parking functions / 
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and g, we say that / <: ^ iff /(v) < g(v) for all vertices v e V. Also recall the parking function 
/ A ^ defined as (/ A ^)(v) := min(/(v), g(v)). Finally, recall dom(/) := {g a parking function : 
8<fl 

Lemma 4.12. Given a non-maximum parking function g, let Fg be the set of maximum parking 
functions that dominate g. Then, g = /\ f^Pg /• 

Proof. It is clear that any non-maximum parking function g is dominated by the meet of all 
maximum parking functions which dominate it. The lemma claims that it is in fact exactly 
equal to the meet. To show this, it suffices to show that for any v, there exists a parking 
function / such that /(v) = ^(v), f{u) > g{u) for all m G V and f{w) > g{w) for exactly one 
vertex w v. Continuing this process, we get a maximum parking function which agrees with 
g on V. Since this is true for all v, we are done. 

To prove the above, note that from the bijection between parking functions and spanning 
trees, we see that any parking function g uniquely corresponds to a mixed orientation of the 
edges of G, where a mixed orientation has some edges directed and others undirected. More- 
over, the directed edges induce a DAG. 

Consider the mixed orientation with respect to g. Let D be the DAG induced by the set of 
directed edges. Suppose there is an edge (w, w) which is not oriented where u and w are distinct 
from v; then in D if there is a directed path from u to w, orient the edge from u to w, else orient 
the edge from w to u. if there is no path either way, orient in any of the two ways. We still have 
a mixed orientation where the directed edges are acyclic. Moreover, this increases the value of 
parking function on exactly one of u or w, and keeps the value the same everywhere else. 

Thus, the non-trivial case is when the only undirected edges are incident on v. Let (m, v) 
be such an edge. Note that we do not want to increase g{v), that is, we want to orient (m, v) 
without increasing v's out-degree. In other words, we want to orient it from u to v. This is a 
problem if and only if there is a path from v to m in D. If so, consider the longest such a path 
P = (v, w, . . . , m), where w is the neighbor of v on P. 

Note that apart from the edge (v, w), there is no path from v to w. If there were such a path, 
it can not use vertices from P, for w "dominates" (in the DAG) all vertices in P other than v. 
Thus, the path from v to w must use "new" vertices making P longer, and thus providing a 
contradiction. 
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Therefore, we can flip (v, w) to (w, v) and not create any cycles. That is, the operation to 
get / is to flip (v, w) to (w, v) and orient (v, u) from v to u. This only increases g(w) but keeps 
everything else the same. □ 



5. Product Graphs and 2„-Parking Functions 

Given two graphs Gi and G2, there is a standard notion of the Cartesian product G10G2 of 
the two graphs. Given a Gi-parking function and a G2-parking function, we define below a 
GiDGi-parking function in a natural way that is symmetric in Gi and G2. 

Definition 5.1. Given G[ = (Vi,E[) and G2 = (V2, E2), the Cartesian product graph G1DG2 = 
(V, E) is defined on the vertex set V = V\ X V2, using the edge set E = E' U E", where 
E' = {{(mi,v),(m2,v)} : [uuU2\ e Eiandv e V2}, and E" = {{(m,vi),(m,V2)} : {vi,V2} e 
E2 and u £ Vi] . 

It is easy to see from the definition that the number of vertices in GinG2 is \Vi\ \V2\, and that 
the number of edges is |V(G2)| |£(Gi)| + |V(Gi)| \E(G2)\. 

5.1. Parking functions on Product graphs. While it seems hard to characterize parking 
functions on GinG2, by simply knowing those on Gi and on G2, the following result gives 
an explicit way to construct a particular family of maximum parking functions on product 
graphs. 

Proposition 5.1. Let fi e 'P(Gi;qi) and f2 e P{G2,q2)- Then f\Uf2 e ;P(GinG2; (^1,^2)). 
where f\Uf2{u, v) = fi(u) + /2(v) + 1, whenever u € Gi and v G G2. 

Further, if fi and f2 are maximum parking functions, then /in/2 is a maximum parking 
function. 

Proof. Let / = /[□/2, G = GinG2, and q = (^1,^2)- To show that / e 'P(G,q), once again 
we make crucial use of Dhar's marking algorithm; in particular, we will make use of the order 
in which the vertices of Gi (and G2) are marked in verifying that f e V{Gi,q\) (and /2 G 
'P(G2,q2), respectively). Using these in turn, we order the vertices in GinG2: let (u',v') <□ 
(m, v), if u' is marked before m in Gi, or if m = u' and v' is marked before v in G2. We now 
prove that / is a valid parking function, by showing that the vertices of GinG2 can be marked, 
using Dhar, in precisely the order given by <□. We begin by noting that f{q) = f(qi,q2) = 
/i(?i) + fiiqi) + 1 = -1. Now consider the vertices of G inductively, using the order given by 
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<□■ By the time the vertex (u, v) is considered, observe that there are at least f(u) + f(v) + 2 
neighbors of (a, v) that have already been marked, since they precede (w, v) in <□: indeed, at 
least f(u) + 1 neighbors of the form {u', v) have been marked, (since the graph induced by fixing 
the second coordinate v is simply Gi), and similarly at least /(v)+ 1 neighbors of the form (u, v') 
have also been marked. Thus (u,v) can be assigned the value f(u,v) = f(u) + /(v) + 1, and 
hence / is a G-parking function. 

Now, suppose that /i and /a are maximum Gi- and G2-parking functions respectively. Then, 
we must show that / achieves the genus of g(G) minus one. This is easy to verify: 

||/in/2l| = |y(G2)|||/i|| + |y(Gi)|||/2|| + |y(Gi)||V(G2)| 

= |y(G2)i |£(Gi)i + |y(Gi)i |£(G2)| - |y(Gi)i |y(G2)i 

= |£(GinG2)| - |y(GinG2)| = g(GinG2) - 1. 

5.2. 2,,-Parking Functions. A quintessential product graph is the n-dimensional cube Q„, 
obtained by taking the product of an edge Qi with itself n times. For the purpose of this 
section, for integer n > 1, we also view the n-cube Q„ as the graph on 2" vertices, which may 
conveniently be labeled by the 2" binary vectors of length n, and with edges between vertices 
whose Hamming distance is one. We are interested in understanding the parking functions on 
Qn with respect to the vertex <5r = (0, 0, . . . , 0). 

Definition 5.2. For n = I, the unique parking function f = on Q\ is canonical. For n> \, 
the parking function f" = /□/□ • • • □/, obtained using the product graph construction, is 

n times 

defined as the canonical Q„-parking function. Further, if a Qn-parking function g is such that 
there exists a bijection (p '■ y(2n)— >y(2n) such that /(v) = g{(p{v))for all v e V{Qn), then we 
say that g is semi-canonical. 



Note that by Proposition l5.ll the canonical 2„-parking function is a maximum parking func- 
tion, and hence the semi-canonical one is also a maximum. 

Example 5.1. Not all maximal Qn-parking functions are semi-canonical. For instance, con- 
sider the Qi-parking function f such that f (000) = -1,/(001) = 1,/(010) = 0,/(100) = 
2, /(Oil) = 0,/(101) = 0,/(110) = 0,/(lll) = 2. It is easy to verify that f is a Q^-parking 
function, but f is not semi-canonical. Since, \\f\\ = 4 = l^l - |y|, / is maximum (and thus 
maximal). 
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Recall the partial order on parking functions, and the notion of dom(/) = {g & PiG, q) : g < 
/}, for a parking function / e !P(G, q). 

Theorem 5.2. Let /" denote the canonical Q„-parking function. Then /"(v) = wgt(v) - 1, 
where wgt(v) is the Hamming weight ( the number of 1 's in the binary representation ) of the 
vertex v e V{Qn). Consequently, if f is semi-canonical, then 

n 

|dom(/)| = ]^^©. 

k=2 

Proof. The proof is by an easy induction on n e The base case is clear, since takes the 
values -1,0. For the induction step, for n > 1, write f"^^ = fnf^, which by Proposition 15. II 
implies that, for u e ViQn) and v e ViQO = {0, 1}, 

/"+i(i<,v) = /» + /i(v) + l. 

Thus for < k < n, the value k can be obtained in /"^^ either by taking a k from /" and adding 
a zero to it (through, /'(O) + 1 = -1 + 1), or by taking ak - I from /" and adding a one to it 
(through, /'(I) = 0+1). By the induction hypothesis, the number of k's is 

completing the induction step. Also note that, by definition, a semi-canonical parking function 
also has the same distribution of integers. 

To complete the proof of the theorem, recall by Proposition l2.2[ that for any v e V(Qn)\{q}, 
if f(v) = k for k e Z>i, then for g < f, we may have ^(v) = 0, 1, . . this implies that there 
are k + I possible values for such a g with respect to v. Since < k < n and each value in the 
range of / is independent with respect to obtaining a dominated function g, we have that 

n 

|dom(/)| = ]~[^G). 

k=2 

■ 

Remark 5.1. Note that it can directly be seen that f is a Q^-parking function, when f is 
defined as /(v) = wgt(v) - 1. Using Dhar's algorithm, starting with q = (0,0, . . . , 0), we may 
proceed by marking vertices level by level ( in the Boolean lattice ): Every vertex v in level k has 
precisely k down-neighbors in level k — \, and they have all been marked, so v can be marked 
and given value k-\. Thus there are precisely {^^ vertices which obtain the value k-\, in such 
a canonical parking function on 2„. 
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Example 5.2. Note that not every semi-canonical parking function is canonical - in satisfying 
the level by level property described in Remark \5.1\ For instance, consider the QT,-parking 
function f such that f {000) = -1, /(OlO) = /(1 00) = /(lOl) = 0, /(OOl) = /(Oil) = 
/(111) = 1, and f (WO) = 2. It is easy to check that f is a Q^-parking function, semi-canonical, 
but not canonical - since, for example, the value 2 is adjacent to a 0. 

We hope the above remarks and examples indicate the difficulty in understanding the max- 
imum parking functions on even a highly structured, symmetric graph such as the n-cube. As 
far as we know, the number of maximum parking functions of Q„ is known only for n < 4. For 
n = 2,3, and 4, this number is 3, 133, and 3040575, respectively. 

6. Diffuse states and acyclic orientations 

As mentioned in the introduction, in the context of chip-firing, the following notion was 
introduced in [1 251 ■ 

Definition 6.1. Given a connected graph G, a function s : G ^ Z'^ is called a diffuse state 
if for every induced subgraph G[A] = {A,E{A)) c G, there exists some u e A such that 
deg(M)U < s{u). Further, let \\s\\ := 2v ■^W- 

Note the (complementary) similarity with the definition of a parking function, by observing 
that deg(M)U < s(u) is equivalent to deg(M) - s(u) < deg(u)\A^. However, we have no special 
vertex such as q. For chip-firing purposes, s{v) may be thought of as the number of chips on v, 
thus \\s\\ denotes the total number of chips in the graph. 

First consider the following algorithm which constructs an acyclic orientation, by using a 
given diffuse state s with \E(G)\ chips. 

Step 1. Since s is diffuse, we may find a vertex v such that deg(v)|G = deg(v) < s{v). 

Step 2. Orient all of the edges incident to v outward; delete v and its incident edges. 

Step 3. The resulting graph is diffuse since it is a subgraph of G; so we may repeat Steps 1-2 
until all edges of the graph are oriented. 

Note that this process gives an acyclic orientation since we cannot orient edges into a vertex 
which has out edges since this vertex has been deleted from the graph. 

Lemma 6.1. If s is a diffuse state on graph G with \\s\\ = \E(G)\, then there is a vertex v with 
s(v) = 0. For every diffuse state s, the set {v : s(v) = 0} is an independent set in G. 
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Proof. The above algorithm which repeatedly removes vertices, removes at least as many chips 
as the edges at each step. Before the last vertex, all edges (hence all \E{G)\ chips) must have 
been removed, which means that the last vertex can has zero chips. For the second part, if 
s{u) = s(v) = 0, and w, v 6 E{G) then the set A = {u, v} violates the diffuse property. I 

With a similar proof, it can also be shown that there exists a vertex w such that s{w) = deg(w), 
under the hypothesis of the above lemma. 

Theorem 6.2. There is a bijection between the set D(G) of diffuse states with \E{G)\ chips and 
the set AO(G) of acyclic orientations of a connected graph G. 

Proof. Let E{G) = m. The proof is based on two injections between the sets, going in each 
direction. First, given an acyclic orientation O e AO(G), define the nonnegative function 
s = sohy letting ^•(v) be the out-degree of v in the orientation O. Clearly, the mapping is one- 
one, since the out-degree sequence uniquely determines an acyclic orientation; also ||5|| = m, 
since the sum of out-degrees equals the number of edges. To see that s is diff"use, simply 
observe that, each induced subgraph G[A] Q G has a (local) source y eA when restricted to the 
acyclic orientation induced on A; such a source y satisfies s(y) > deg(j)U, since the out-degree 
of y is at least the degree deg(y)U. 

For an injection in the other direction, we make use Lemma [611 Given a diff'use state s with 
m chips, we construct an acyclic orientation, by constructing a DAG: Lemma [6T| guarantees 
the existence of sink(s); so we construct the orientation, by (i) repeatedly removing the current 
set of sinks, and (ii) subtracting a chip from each in-neighbor of a removed sink. It is easy 
to see that the updated function s at each step is still a diffuse state on the remaining graph. 
Note that this construction is one-one: for s, s' different diffuse states, simply consider the first 
time the current sets of sinks (in the DAGs) differ, when we start with s versus s'; since the 
underlying graph is the same, there must be such a time whenever s s' . ■ 

Remark 6.1. Thanks again to an anonymous referee of^, we learnt the following: that diffuse 
states are also in 1-1 correspondence with the critical configurations of the chip-firing game of 
Bjdrner-Lovdsz-Shor [fTT| . and that the first part of the above proof appears as Theorem 3.3, 
part (b) in [[TT|. In addition, it is easy to go between the above theorem and Theorem \3.1[ 
simply add a new vertex q adjacent to every vertex in G and obtain a new graph G; then 
G-parking functions (with respect to q, say) correspond to diffuse states in G, and acyclic 
orientations with a unique source at q in G correspond to acyclic orientations ofG. 
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7. Concluding Remarks 

Soon after the completion of this work, Igor Pak kindly pointed us to the work of Olivier 
Bernardi [|3, where bijective proofs are derived for interpretations of each of the evaluations 
of the Tutte polynomial TG{x,y), fox < x,y < 2, in terms of orientations. A key to this seems 
to be a nice combinatorial embedding of (the edges of) the graph and a rewriting of the Tutte 
polynomial using notions of internal and external embedding activity. 

In addition to the questions mentioned in the previous sections, several challenging problems 
remain open. Given an arbitrary graph G, it is a classical open problem in the topic of Markov 
chain Monte Carlo (MCMC) algorithms ll26l . to efficiently generate an acyclic orientation 
uniformly at random from the set of all such orientations. Due to the observations above, a 
closely related problem would be to generate at random an acyclic orientation with a uniquely 
identified sink (or source), or equivalently, to sample uniformly from the set of safe spanning 
trees of G. The MCMC technique suggests the following natural approach to this problem: 
it is well known (see for example, IfT^ [T9]| ) that the so-called bases exchange walk provides 
an efficient way to sample uniformly from the set of all spanning trees of a given graph G. 
However it remains to be seen whether (and how) restricting such a random walk to the set 
of safe trees affects the mixing time of the walk - the time by which the walk converges to its 
steady state distribution, uniform on the set of safe trees. Given that the exact enumeration of 
the number of safe trees of 2„ is also open, an interesting first step might be to analyze such a 
walk on the trees of 2„. 

Independent of the above approach, other ways of providing asymptotically accurate esti- 
mates, for large n, of the number (or even the logarithm of the number) of maximum parking 
functions on Qn , remains interesting and presumably a challenging exercise. 
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